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Tree Calculus for Bivariable Difference Equations 
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psj _ Abstract. Following Poupard's study of strictly ordered binary trees 

with respect to two parameters, namely, "end of minimal chain" and 

rS ,' "parent of maximum leaf" a true Tree Calculus is being developed to 

.^ I solve a partial difference equation system and then make a joint study 

. of those two statistics. Their joint distribution is shown to be symmetric 

0^ ' and to be expressed in the form of an explicit three-variable generating 

function. 

o' 

^ ' 1. Introduction 

r^ ' The triangle of numbers 

-(— > I 

^ ' /o(l) 1 

/i(l) /i(2) /i(3) 1 

/= /2(1) /2(2) /2(3) /2(4) /2(5) = 12 10 

/3(1) /3(2) /3(3) /a (4) /a (5) /a (6) /a (7) 4 8 10 8 4 

> : /4(1) /4(2) /4(3) /4(4) /4(5) /4(6) /4(7) /4(8) /4(9) 34 68 94 104 94 68 34 

^ : Table 1.1 

■^ ! appears in Sloane's On-Line Encyclopedia of Integer Sequences [S106] 

under reference A008301 and is called Poupard's triangle. As shown by 
Christiane Poupard [Po89], / = (/„(to)) (n > 0, 1 < m < 2n + 1) is the 



Cf) • unique solution of the finite- difference equation system 

(1.1) A2/r.M + 2/n-i(m)=0 (n>l, l<m<2n-l), 

is^ . where A stands for the classical finite- difference operator (see, e.g., [Jo39]) 

^: (1-2) A/„(m):=/^(m + l)-/„(m), 

so that 

(1.3) A^Uim) = Uim + 2) - 2/^(m + 1) + /„(m), 

when taking /o(l) = 1; ^(1) = 0, /„(2) = Y, /n-i("i) for n > 1 as initial 
values. Let „ , "^ 

(1.4) tan« = X:^^--^T2._i 

~ S 'i 7 Q 

U U U U U 
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be the Taylor expansion of tanw, the coefficients T2n+i {n > 0) being 
caUed the tangent numbers (see, e.g., [Ni23, p. 177-178], [Co74, p. 258- 
259]); Poupard further shows that each row sum /„(•) := /n(l) + /n(2) + 
■ ■ ■ + /n(2n-|-l) is equal to the integer T2n+i/ '2"' (n > 0), that is, reporting 
to Table 1.1: 1, 1, 4, 34, 496,. . . 

Finally, on the set l2n+i of strictly ordered binary trees with (2n + 1) 
vertices (see Definition 1.3), she defines two statistics "eoc" ("end of 
minimal chain" ) and "pom" ( "parent of the maximum leaf" ) , to show 
that both statistics "eoc" and "pom-|-l" are equally distributed on each 
set T2n+i, and furthermore, 

(1.5) #{t G T2n+i:eoc(t) = k + l} = #{t G T2n+i:pom(t) = k} = U{k) 

for all k; in particular, ^%2n+i = T2n+i/'^"'- 

The purpose of this paper is to calculate the joint distribution of the pair 
(eoc, pom) on each set l2n+i and to derive its properties, in particular its 
symmetry. To achieve this, we first introduce a sequence (M^ = {fn{m^ k)) 
of (2n) X (2n)-matrices (n > 1) with nonnegative integral entries, called a 
Delta sequence, defined by a system of finite difference equations, verifying 
certain initial conditions. Then, we show that each entry /n("^, k) is equal 
to the number of trees t from %2n+i such that eoc(t) = m and pom(t) = k. 

It is convenient to consider the following four triangles of each square 

{(m, k) : 1 < m, k < 2n}: 

L^J:'^ :={2<k + l<m<2n- 2}; L^^^ := {4 < k + 3 < m < 2n}; 
UJ^-'' :={2<m + l<k<2n- 2}; ui^^ := {A < m + 3 < k < 2n}. 

By convention, /n(?Ti, k) :=0 if (m, k) ^ [1, 2n] x [1, 2n]. The partial differ- 
ence operators A, A, act as follows on the entries of the matrices M„: 

m k 

(1.6) Afn{m,k) := /^(m + l,A;) - fn{m,k); 

m 

(1.7) A /,(m, k) := fn{m, k + l)- f^{m, k). 

k 

They serve to define the recurrence relations: 

(Rl) AVn(m,A;) + 2/^_i(m,A;) = ((m,A;)GLW); 



m 



{R2) AVn(m,A;) + 2/^_i(m,A;) = ((m. A;) G t/^). 

k 

Finally, denote the row and column sums of M^ = {fn{T^i k)) by 



fr,{m,.):= 2_^ fn{m,k); (1 < m < 2n); 

l<k<2n 

fn{'.k):= V fn{m,k) {l<k<2n). 



l<m<2n 
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Definition 1.1. A sequence of matrices (M„) {n > 1), where each 
matrix M„ = {fniin, k)) {1 < k,m < 2n) has nonnegative integral entries, 

having only O's along its diagonal, and such that Mi '•= [ ^ n ) ' ^^ ^^^^ ^° 

be a Delta Sequence, if for n > 2 both recurrence relations (-R1) and (-R2) 
hold, together with the initial conditions: 

[1 1) for n > 2 the (2n)-th column, Col2n, of M^ is the zero-column; its 
the (2n — l)-st column, Col2n-i, is equal to 

/„_i(l,.), /n-i(2,.), •••, /n-i(2n-2,.), 0, 0, 
when read from top to bottom. 

(/ 2) the (2n)-th row of Mn is also equal to 

/„_i(l,.), /n-i(2,.), •••, /n-i(2n-2,.), 0, 0, 
when read from left to right; its {2n — l)-st row is equal to: 

/n-l(l, •) + /n-l(», 1), /n-l(2, .) + /n-l(»,2), ... , 

/,_i(2n-3,.)+/n-i(.,2?i-3), /,_i(2?i-2,.) + /n-i(.,2?i-2), 0, 0, 

when read from left to right. 

In the above definition the entries of the matrix M„ are derived 
from Mn-i by first applying rules (/I) and (/2) and letting the diagonal 
be null; then, starting from m = 1 up to m = 2?i — 3, for each k 
from 2n — 3 down to m + 1, evaluate fn{rn,k) with equation (-R2): 
/n(TO, k) —2fn{m, k + l) + fn{m, k + 2) +2fn-i{m, k) = 0, the coefficients 
fn{m,k + 1), fn{m,k + 2) and fn-i{m,k) being already calculated. 
Exchanging the roles of m and k the upper entries are obtained by using 
equation (-R1). Accordingly, (-R1), (-R2), (/I), (/2) uniquely determine 
the Delta Sequence (M„ (n > 1). 



Calculation of the first matrices. 


First, /i 


that 

M2 = 


? 
1 

Vo 


? 0\ 

1 

1 

1 0/ 



0, /i(2. 



1, so 



by rules (/I) and (12). The remaining entries are obtained by rule (-R2): 
/2(1,2) -2/2(1,3) + /2(1,4) + 2/1(1,2) = /2(l,2)-2x + + 2x0 = 0, 
so that /2(1, 2) = 0; then, by rule (Rl): ^(2, 1) - 2/2(3, 1) + /2(4, 1) + 
2/1(2, 1) = /2(2, 1) -2x1 + + 2x1 = 0, so that /2(2, 1) = 0. Thus, 



Mo 
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M, 



^0 ON^ 
12 10 
110 4 2 
2 3 4 10 
13 3 10 

\0 1 2 1 O/ 



M, 



/OOOOOOOOX 

0048 10 840 

4 4 16 20 16 8 

8 12 16 28 20 10 

10 18 24 28 16 8 

8 18 24 24 16 4 

4 12 18 18 12 4 

\0 4 8 10 8 4 0/ 



M. 



/ 



34 
68 
94 



000000000 

34 68 94 104 94 68 34 

34 136 188 208 188 136 68 

102 136 274 296 262 188 94 

162 222 274 352 296 208 104 



104 198 276 330 352 274 188 

94 198 282 330 330 274 136 

68 162 240 282 276 222 136 

34 102 162 198 198 162 102 34 

\ 34 68 94 104 94 68 34 



94 

68 








34 

/ 



Fig. 1.2: the first matrices M„ 



and/2(l,.)=0, /2(2,.) 
are displayed in Fig. 1.2. 



I7 /2(3, •) = 2, 72(4, .) = 1. The next matrices 



The previous definition of a Delta Sequence, based on the two relations 
(i?l), (-R2) and the two initial condtions (/I), (-^^2), can be symbolized 
by the square on the left in Fig. 1.3, as relation (-R1) (resp. (-R2)) acts on 
the entries of the lower (resp. upper) entries of the matrix M^, and initial 
conditions (/I) and (/2) refer to the last two columns Col2n-i, Col2n and 
rows Row2n-i, Row2n of M^, respectively. Other initial conditions will be 
stated in Section 8. We just mention a second one, materialized by the 
square on the right in Fig. 1.3. 



Col2„_l Col2„ 



RoW2„_l 

Row2„ 




Coll C0I2 
Rowl i^' 
Row 2 



>(i?2) (i?4) 




(Rl) 



Fig. 1.3: Definitions 1.1 and 1.2 
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Definition 1.2. A sequence of matrices (M„) {n > 1), where each 
matrix M„ = {fni^Ti, k)) {1 < k,m < 2n) has nonnegative integral entries, 

having only O's along its diagonal, and such that Mi '•= [ ^ n ) ' ^^ ^^^^ ^° 

be a Gamma Sequence, if for n > 2 both recurrence relations 

(i?3) A''U{m,k) + 2U-i{m,k-2) = {{k,m) e U^^^); 



m 



{R4) A''U{m,k) + 2U-i{m-2,k) = ((A;,m) G 4^)); 

k 

hold, together with the initial conditions: 

(/3) for n > 2 the first row, Rowi, is the zero-row; the second row, 
R0W2, is equal to 

0, U-i{l, .)(= 0), /^_i(2, .),..., fu-i{2n - 2, .),fn-i{2n - 1, .)(= 0); 

when read from left to right. 

(14) the first column, Coli, of M^ is also equal to 

0, /^_i(l, .)(= 0), /^_i(2, .),..., /n-i(2n - 2, .), /n-i(2n - 1, .)(= 0) 

when read from top to bottom; the second column, C0I2, is equal to 

0,0,/„_l(2,.)+/n-l(l,.), /n-l(3,.)+/n-l(2,.), ••• 

fn-i{2n - 2, .) + /n-i(2n - 3, .), /n-i(2n - 2, .), 
when read from left to right. 

Using the same reasoning as for Definition 1.1 it is seen that the Gamma 
Sequence is uniquely defined. The fact that Delta and Gamma Sequences 
are identical will be a consequence of the further theorems (c/. Section 6). 
Next, comes the combinatorial set-up on which all calculations will be 
made. 

Definition 1.3. An n-laheled strictly ordered binary tree is defined by 
the following axioms: 

(1) it is a labeledtree with n nodes, labeled 1, 2, . . . , n; the node labeled 1 
is called the root; 

(2) each node has no child (it is then called a leaf), or two children, 
their order being immaterial (it is then called an interior node); 

(3) when getting along each path from the root to each node, the node 
labels are in increasing order. 

Each strictly ordered binary tree has an odd number of vertices, say, 
2n + 1, with n interior nodes and n + 1 leaves. Let ^2n+i denote the set 
of all strictly ordered binary trees with {2n + 1) nodes. When giving an 
orientation (left or right) to each child of each of the n interior nodes in 
a tree t G T2n+i, we generate 2"- planar strictly ordered binary trees (also 
called "arbres binaires croissants complets" by Viennot [Vi88, chap. 3, 
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p. 111]). It is known that the latter are equidistributed with the alternating 
permutations of order (2n + l), so that their number is equal to the tangent 
number T2n+i, a result that goes back to Desire Andre [Anl879, Anl881]. 
Accordingly, 



(1.8) 



#2:^ 



2n+l 



T, 



2n+l 



/2'". 



Let t G T2n+i (n > 1). If a node labeled a has two children labeled b 
and c, define niina := min{6, c}; if it has one child 6, let mina := b. The 
minimal chain of t is defined to be the sequence ai —)■ a2 —)■ as —)■■■■ —)■ 
aj-i — )■ Oj, with the following properties: (i) ai = 1 is the label of the 



root; (ii) for i = 1, 2, . . 
interior node and ai^i 



,j — 1 the (z + l)-st term ai_|_i is the label of an 
: mintti; (iii) aj is the node of a leaf. Define the 
"end of the minimal chain" of t to be eoc(t) := aj. If the leaf with the 
maximum label (2n + 1) is incident to a node labeled k, define its "parent 
of the maximum leaf" to be pom(t) := k. 





Fig. 1.4 

The minimal chain of the tree t displayed in Fig. 1.4 on the left is 
1 — 7- 2 — 7- 3 — 7- 7, so that eoc(i) = 7 and the parent of its maximum 
leaf (equal to 2n + 1 = 9) is pom(t) = 4. Also, the parent of the 
maximum leaf in the tree t' on the right is equal to pom(t') = 6. To 
go from t to t' replace the labels 1,2,3,7 of the minimal chain by the 
labels 2 — 1,3 — 1,7 — 1,9 = 1,2,6,9 and change each other label a 
by a — 1. We then have eoc(t) = pom(t') + 1 = 7. This illustrates 
the construction of the bijection t i— ?> t' of %2n+i onto itself with the 
property eoc(t) = pom(t') + 1, which was described in [Hal2]. Thus, 
the statistics "eoc" and "pom+1" are equidistributed on each %2n+i, as 
already mentioned in (1.5), a property first proved by Poupard [Po89] by 
analytic methods. The main results of this paper are stated next. 

Theorem 1.1. Let (M„ = {fn{rn, k)) {n > 1) be the Delta sequence, as 
introduced in DeSnition 1.1. Then, for aU n > 1 and 1 < m,k < 2n 

(1.9) #{t G l2n+i : eoc(t) = m, pom(t) = k} = fn{m, k). 
A first consequence of Theorem 1.1 and (1.8) is the identity: 

(1.10) 5^/„(m,A;)=T2,+i/2" (n > 0). 



6 
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Theorem 1.2. Let (M„ = {f^{m,k)) (1 < m, /c < 2n)) {n > 1) be the 
Delta sequence. Then, the matrices M„ are symmetric with respect to 
their counter-diagonals: 

(1.11) fn{m,k) = f^{2n + l-k,2n + l-m,) (1 < /c, m < 2n). 

The proofs of those two theorems and also of further properties will 
be based on the geometric properties of the strictly ordered binary trees. 
To this end, we adopt the following notation and convention: for each 
triple (n,m, /c) let l2n+i,m,k (resp. l2n+i,m,,, resp. T2n+i,.,fc) denote the 
subset of T2n+i of all trees t such that eoc(t) = m and pom(t) = k (resp. 
eoc(t) = m, resp. pom(t) = k). By convention, designate those families and 
their cardinalities by the same symbol and also the matrix of the integers 
^2n+i,m,k by Mat(T2n+i)- Our plan of action will be to show that the 
sequence (Mat(T2n+i)) {n > 1) is identical to the Delta Sequence. 

In Sections 2-5 it will be shown that, when replacing each fniin, k) by 
'^2n+i,m,k the initial conditions (/I) and (/2), the two finite- difference 
equations systems (-R1), (-R3), the two finite-difference equations systems 
{R2), {R4) and the initial conditions (/3) and (/4)) hold. This will 
complete the proofs of Theorems 1.1 and 1.2, as done in Section 6. Further 
properties of the matrices Mat(T2n+i) (and then matrices M^) will be 
given in Section 7. Finally, several other equivalent definitions of the Delta 
sequence will be mentioned in Section 8. In Section 9 we conclude the paper 
by calculating the generating functions for the /n(?Ti, k) in the following 
forms. 

Theorem 1.3. The triple-exponential generating function for the lower 
triangles of the matrices M^ is given by 

rn — k — l k — 1 2n — m 

(1.12) Y: ^-("^-^) (^-fc-i)!(fc-l)! (2^ 

cos(v^a^) +cos{y/2y) cos(\/2; 



to)! 



^ ^fx + y + z 
2 cos^ 1= — 

V ^/2 

Theorem 1.4. The triple exponential generating function for the upper 
triangles of the matrices M^ is given by 

™2n— fc „ k — m—1 ^m — 1 



(1.13) V U{m,k)— — — — - 

^-^ 2n — fc ! A; — m — 1 ! m — 1 



2<m+l<k<2n-l 

= sin(v2a;) sin(v2; 



)! 
1 



2fx + y + z 
2 cos ' 



V2 
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2. The initial conditions (/I) and (/2) 

In this section and the next ones we make the convention that whenever 
a leaf is deleted from a tree, the edge linking the leaf to the tree is also 
deleted. 

For verifying that the matrices Mat(T2n+i) have only zero in their 
diagonals, it suffices to show that 'Z2n+i,m,m = (or is equal to with our 
convention). This is true, because if t G l2n+i and eoc(t) = pom(t) = m, 
the node (2n + 1) has a parent equal to m. Consequently, m cannot be 
the end of a minimal chain. Hence, the previous subset is empty. 

Theorem 2.1. The initial conditions {1 1) and (12) hold for the matrices 
Mat(T2n+i) when fn{m^*) and fn{*,k) are replaced by l2n+i,m,. a.nd 
'^2n+i,,,k, respectively. 

Proof. (/I) First, the (2n)-th column of the matrix Mat(T2n+i) has 
zero entries only, as pom < 2n — 1. Next, each tree from '^2n+i,m,2n-i 
(1 < TO < 2n — 2) must contain the subtree 

2n, 2n+l 



Hence, ^2n+i,2n-i,2n-i is empty, for (2n — 1), being an interior node, 
cannot be the end of the minimal chain. Also, ^2n+i,2n,2n-i is empty, 
for the sibling of (2?i — 1) is neccessarily less than (2n — 1), so that the 
minimal chain cannot go through (2?i — 1) and reach (2n). Furthermore, 

'J2n+l,l,2n-l = 2^2n-l,l,. = 0, aS eOC > 2. 

In the remaining cases, that is, 2 < to < 2n — 2, removing the two leaves 
(2n), (2n + 1) transforms each tree from ^2n+i,m,2n-i onto a tree from 
^2n-i,m,, in a bijective manner. Such a transformation may be illustrated 
by the diagram: 




2n+l 
2n-l m^r 7 2n-l 




1-^ 
1 1 

Hence, the (2n — l)-st column of the matrix Mat(T2n+i) reads: 

(2.1) "3^2^-1,1,. 5 2^271-1, 2,.7 ••• 7 2^271-1,2^-2,. 7 0, 0. 

from top to bottom. 

(/2) For the (2n)-th row of the matrix Mat(T2n+i) note that 
2^2n+i,2n,i = for n > 2. When k > 2 each tree from l2n+i,2n,k must 
contain the subtree 2n« 2n+i 
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By (/ 1) we then have: l2n+i,2n,2n-i = '^2n+i,2n,2n = 0. For the remaining 
cases 2 < k <2n — 2 we can set up a bijection of%2n+i,2n,k onto T2n-i,A;,. 
by removing the two leaves 2n and (2n + 1), as illustrated by the next 
diagram. 

2n+l 





1-^ 
1 1 

Note that the node k becomes the end of the minimal chain. Thus, the 
(2n)-th row of the matrix %2n+i is also equal to (2.1) read from left to 
right. 

Finally, consider the {2n — l)-st row of %2n+i- In an obvious manner, 

l2n+l,2n-l,2n-l = l2n+l,2n-l,2n = . When 1 < /c < 2?1 - 2, the 

trees from the sets l2n+i,2n-i,k fall into two categories '^in+i 2n-i k ^^^ 
^2n+i 2n-i if ^^ ^^^ ^^^^ Category the trees contain the subtree 

2n+l • 2n-l 2n • 2n-l 

^^ : in the second one, the subtree ^^ [a ^ k). 

First, note that l2^_^i 2n-i i = '^2n-i,i,» = 0. When 2 < k < 2n — 2, 
removing the two leaves (2n + 1), (2n — 1) and replacing the node label 
(2n) by (2n — 1) maps 'Z2n+i 2n-i k onto T^2n-i,k,» in a bijective manner. 

When 1 < k < 2n — 2, removing the two leaves (2n), (2n — 1), and 
replacing the node label (2n + 1) by {2n — 1) maps ^2n+i 2n-i k onto 
'^2n-i,,,k in a bijective manner. Thus, the {2n — l)-st row of l2n+i reads 

'^2n-l,l,, +^2n-l,,,l, '3^2n-l,2,. + 2^2n-l,.,2, ■■• , 

2^271-1,2^-3,. + 2^2n-l,.,2n-37 '2^2n-l,2n-2,. + 2^2n-l,.,2n-2, 0, 0. [] 

3. Tree Calculus for the relations (-R1) and (-R3) 

In the following Tree Calculus subtrees (possibly leaves) are indicated 
by the symbols "O?" "V"? or "O." The end of the minimal chain in each 
tree is represented by a bullet "•." Letters occurring below or next to 
subtrees are labels of their roots. For instance, the symbols 





designate the families of all trees t from the underlying set %2n+i having 
a node labeled h [in short, a node 6], parent of both a subtree of root a 
and the leaf m, which is also the end of the minimal chain; moreover, the 
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symbol on the right has the further property that the node labeled c does 
not belong, either to the subtree of root 6, or to the path going from root 1 
to b. In the sequel, the letter "m" is always used to designate the end of 
the minimal chain, unless explicitly indicated by a letter next to •. 

Our Tree Calculus consists of two steps: (a) decomposing the sets 
'^2n+i,m,k into Smaller subsets by considering the mutual positions of the 
nodes to, (to + I), {m+2) (resp. k, (/c + 1), (/c+2)); (b) setting up bijections 
between those subsets by a simple display of certain subtrees, as done in 
(3.1). 

For instance, 





(3.1) Cs := I and Di := 

may be regarded as two subsets of %2n+i- To each pair (^^2'*^)' ^^^ 
nodes of "O" being all greater than or equal to {m + 2), there correspond 
a unique tree from C3 and a unique tree from Di. This clearly defines a 
bijection of C3 onto Di. 

Those two principles (a) and (b) will be applied in the proofs of the 
next two theorems 3.1 and 4.1. 

Theorem 3.1. If (m, k) belongs to L^n^ U U^^^ = {2<k + l<m< 
2n-2}U{4<m + 3</c< 2n}, then 

m+l\/rn+2 

(3.2) A '^2n+l,m,k + 2 I =0, 

m 

with the understanding that the second term on the left-hand side repre- 
sents twice the set of all trees from %2n+i,m+i,k with the further property 
that m is the parent of both (m + 1) and (m + 2). 

Proof. The decomposition 

^2n+l,m,k = I + 

means that in each tree from T2n+i,m,A; the node {m + 1) is, or is not, the 
sibling of the leaf m. In the next decomposition the node m is, or is not, 
the parent of the leaf (to + 1): 






m+l 
'^2n+l,m+l,k — I + I {b ^ m) . 



10 
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Under the transposition (m,?B + l) the node labeled k remains unaffected, 
because /c < m — 1 if (m, k) E Ln and m + 3 < /c if (to, k) E Un \ so 
that the parent of (2n + 1) remains k. Thus, the transposition establishes 
a one-to-one correspondence between the two second terms. Hence, 



A^T 



2n-\-l,Tn,k 



CI: 



2n+l,m+2 



k—% 



2n+l,m+l,k) 



C^T: 



2n+l,m+l 



,k — % 



2n-\-l,m,k) 



m+2 




m-\-l 




o 

m+2 



m+1 





o 

m-\-l 



+ 

:= A - B - C + D . 

Depending on the mutual positions of nodes to, (m + 1) and (to -|- 2) 
the further decompositions prevail, as again, k remaining still attached to 

(2?i + l): 

m+T\/^\ — I m+2^ 
m+l\y^ m+1' 

A= f^ + [ I ,m]:=Ai + A2; 

m+l\ /"1+2 m,+ l\ Z-^+S m+1 

B= r + Y"^ + [ 





o 

'm+2 



to] := Bi + B2 + Bs; 



TO + 2] 




+ [ 



m + 2] ■=Di+D2. 



The permutation ( ^2 ^ m+i ) establishes a one-to-one correspon- 
dence between A2 and C4 (resp. between S3 and 02). On the other hand, 
Ai = B2 + C2 = 2 i?2, because the two subtrees in B2 = C2 are incident 
to the same node, while in Ai they are incident to two different nodes. 
Finally, C3 = Di as explained in (3.1). Making all the proper cancellations 



in the sum A %2n+ 



l,m,k 



A -S-C + L* we get identity (3.2). □ 



Corollary 3.2. The relations (-R1) and (-R3) hold for the matrices 

Mat(T2n+i), that is, 

(3.3) A2 %2n+l,m,k + 2 l2n-l,m,k = 0, if (to, k) E L^^); 
m 

(3.4) A^l2n+l,m,k+2l2n-l,m,k-2 = 0, if (to, /c) E U^^l 
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Proof. For (3.3) change the second term in (3.2) as foUows: remove the 
two leaves (m + 1), (m + 2), and subtract 2 from aU the remaining nodes 
greater than (m + 2): the term becomes l2n-i,m,kj as the node (2n + 1) 
becomes (2n — 1) and is stiU hnked to k. For (3.4) do the same changes, 
but this time, as to + 3 < /c, the edge going from k to (2n + 1) becomes 
and edge going from (k — 2) to {2n — 1). G 

4. Tree Calculus for the relations (-R2) and (-R4) 

Theorem 4.1. If {m,k) belongs to L^n^ U U^^ = {4</c + 3<m< 
2n} U {2 < m + 1 < A; < 2n - 2}, then 



(4.1) 



k 



2n+l 



2n+l,m,fc 



+ 2 [ ? , m] + 2 




= 0, 



with the understanding that the second term on the left-hand side is twice 
the number of all trees from %2n+i,m,k+i with the further property that 
{k + 2) is a leaf incident to {k + 1), itself incident to k, the end m of the 
minimum chain being outside the subtree of root k. 



Proof. First, 



fc+i 




2n+l 



2^2n+l,m,fc = [ I" ,w] + [ r ,m, /C + 11 + 
:= Ai + ^2 + A3 + A4 




2n+l 



+ [ 




2n+l 



/c + 1] 



meaning that each tree from %2n+i,m,k has one of the four forms: either 
/c + 1 is incident to /c, or not, and m is outside or not the subtree of root /c; 
furthermore, the leaf m is the end of the minimal chain. 
Using the same dichotomy. 



^2n+l,m,k+l 



Consider the subsets A'^ : 





2n+l rn 



2n+l \/ [^ ^ ,2n+l 
k+1 

fc+1 7 1, Ya- , r Y^+1 7 n 

,m, /c] + r + [ I , /c] 



2n+l 




+ 1 k+1 

of A4 and S^ := I of B2. The 

transposition (/c, /c + 1) maps A4 \ A'^ onto S4 and A2 onto S2 \ -B2 in a 
bijective manner. 



12 



TREE CALCULUS FOR DIFFERENCE EQUATIONS 



Hence, 

'^2n+l,m,k + l — '^2n+l,m,k = (-^i — Ai) + ((i?2 — B2 — A2) + -B2) 

+ (S3 - As) + (S4 - (A4 - A'^) - A'^) 
= Bi-Ai+B'2+B3-As-A'^ 

.2n+l Q ,2n+l 

,2n+l 





2n+l 



and also 



2n+l 



^2n+l,m,k+2 ~ ^2n+l,m,k+l — [ 




+ 



2n+l 

2n+l X J^ 

fc+2\ / fc+2\ / A;+l\ yfc+2 

fc+l Yfc+1 




:=Di-Ci + I)^+I)3-C3-Cl. 



Thus, 



A T2n+l,m,fc — ('J^2n+l,m,fc+2 " '3^2n+l,m,fc+l) — ['^2n+l,m,k+l ~ ^2n+l,m,k) 

-Bi + Ai-B!,-Bs + As + A'^. 
The further decompositions of the components of the previous sum 
depend on the mutual positions of the nodes /c, {k + 1), (fc + 2). First, 
evaluate the subsum: Si := Di — Ci — Bi + Ai using the decompositions: 



2n+l 



Di = [ 




2n+l 



m\ 




2n+l 



m,k] + [ 




,m 



1^1,1+1^1,2; 



2n+l 




m\ 
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,2n+l C) ,2n+l 



Bi 





2n+l 



[ r , m, A; + 2] + [ 




2n+l 



'^ , m] + [ 




k+2 

m] + [ \k ,771] 

Bi^i + -Bi,2 

+ Bis + Bi 4 + -Bi 5; 




^1,1 +^1,2- 



Also, let 



D'1.1 ■■ 




2n+l 

k+2\y y« 
k+l\ /k 



C'ly. 




The permutation ( ^.^ ^ a;+i) '^^Ps -Di,i \ -Di 1 onto i?i,i and C\^\ \ 
C{ 1 onto Ax^x. Hence, Di^i = Si^ + D^^, Ci,i = Ai^i + C'xx- Moreover, 
DI2 = 2Si,3, Cx,2 = Ax,2. Bx,2 = 'Bx,4.'Bx,3 = Bx,^'. Altogether, 

Sx=Dx-Cx-Bx+Ax = {Bx,x + D'xx + 2Bx,3)-{Ax,x + C'xx + Ax,2)- 
{Bx,i + Bx,2 + Si,3 + Bx,2 + Si,3) + (Ax,x + Ai,2). Thus, 

(4.2) Sx = -2Bx,2 + D[^^-Cl^. 

Next, evaluate the sum ^'2 := D2 + ^3 - C'3 - C4 - S^ - S3 + A3 + A^ 
by decomposing each of its components. We have: 



D' 




-^2,1 + -^2,2J 



2n+l 




1^3,1+^3,2; 
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"1 "I i; ,2n+l 

2n+l KJ ,2n+l 

fc+2\ / fc+2\ / k+V 





C'a.i + C'3,2; 



C^ 




2n+l L^ 2n+l W 2n+l 

n \/ n \/ fc+2 

fc+l\ /fc+2 A;+l\ /fc+2 fc+i\ / fc_|_i 

Tfc 



/c] + 

^4,1 + ^4,2 + ^4,3! 




+ 




B', 




'■= -^2,1 + ^2,2 + ^2,3; 




A; + 2]+ r + 
-83,1 + -B3 2 + -63,3 + -63,4; 




,2n+l IJ .2n+l 

fc+l\ / fc+l\ / 

^3= I' = [ I' , k + 2] + 



k+2 




2n+l 



^3,1+^3,2; 



A' 





k + 2] + 

^4,1 + ^4,2 + ^4,3- 



Within the sum 152 there are numerous canceUations we now describe. 

(a) Components of the form [t,k] or [t, /c + 2], where t is a subtree, 
whose root is labeled. There are four of them: -Da,!, —6*3,1, —-63,1, ^3,1- 
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Consider the subsets: m 

,2n+l 

fc+r 




-S3, 1,1 '■— I 7 ^3,1,1 • = 




k+2 



oi B31 and ^3^1, respectively. The permutation {j^j^2 k k+i) i^^ips Ds,i 
onto Ss,! \ S3, 1,1 and Cs,! onto ^3,1 \ A3, 1,1. Hence, 1)3,1 - C3,i - ^3,1 + 

^3,1 = (-^3,1 — -S3, 1,1) — (^3,1 — ^3,1,1) — Bs,l + As,i = —-63,1,1 + ^3,1,1- 

(b) Components of the form [t,k] or [t,k + 2], where t is a subtree, 
whose root is not labeled. There are four of them: D2 1, — C4 i, — i?2 1, 
A'^-i- Again, the permutation ( ^ ^ ^ ^^^ ) maps D'2^ onto -62,15 and C'^^ 
onto A4 ^. Hence, -D2 1 — -B2 1 = — C4 ^ + A4 ^^ = 0. Their sum vanish. 

(c) Components represented by a tree t, whose root is unlabeled. There 
are four of them: —B22, — -^2 3, — ^4 2' ^4 3- As -B2 2 = ^4 2' ^^^ 
contribution of those components to 5*2 is then —B2 3 + A4 3. 

(d) Components represented by a tree t, whose root is labeled. There 
are nine of them: 1^2,25 ^3,2, -(^3,2, -C42, -(^4,3 --83,2, --83,3, --83,4, 
A3, 2. By simply comparing the subtree contents we have: D22 — C'3,2 = 
-^3,2 + ^3,2 = 0, 1^3,2 - (C;3 + -83,4) = and C'^^2 = 53,3- The 
contribution of those terms is then —2 C'^ 2- 

Hence, 5i + 52 = (-2Si,2 + I)l,i-C(,i) + ((-S3,i,i + A3,i,i) + (-S^,3 + 
^4,3) + (-2^;2))- As D[^^ = S^^3, C[^, = ^3,1,1 and ^3,1,1 = ^^,3, 
get S'l + 5*2 = -2 Si,2 - 2 C4 2, that is, 

A^%2n+i,m,k-2[ P ,m] -2 ? . D 

k 

Corollary 4.2. The relations (-R2) and (-R4) hold for the matrices 
Mat(T2n+i)j t-hat is, 

(4.3) A^ %2n+l,m,k + 2 T2n-l,m-2,fc = 0, if (m, k) E L^f ^ 

k 

(4.4) A^%2n+l,m,k+'2%2n-l,m,k = 0, if {m,k) G t/ « . 

k 

Proof When (m, k) belongs to Ln , the second term in (4.1) is in 



we 






2n-l 

bijection with 2[ I "^ , to — 2] and the third one with 2 I "^ , that is, 

the set of trees in which the end to — 2 of the minimal chain is outside 
(resp. inside) the subtree of root k. The sum of those two terms is then 

2T2n-l,m-2,fc- 
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When (to, k) belongs to t/A , the third term of (4.1) vanishes and the 




2n-l 

second one is in bijection with 2[ 1 ,771], equal to 2T2n-i,m,fc- D 

5. The initial conditions (/3) and (/4) 

Property 5.1. Initial conditions (/3) and (/4) hold for the matrices 

Mat(T2n+i). 

Proof. (/ 3) The first row of each matrix Mat(T2n+i) is obviously the 
zero-row, as 1 can never be the end of the minimal chain. For the second 
row note that for n > 2 the set %2n+i,2,k is empty when /c = 2 or 2n. Also 
the set %2n+i,2,i is empty, for 2 and (2n + 1) can be both children of the 
root only when n = 1. 

Let 3 < /c < 2?i — 1. As illustrated by the diagram 

2n+l 





1 

each tree from T2n+i,2,A: is transformed into a tree from '^2n-i,,,k-2 by 
deleting the two nodes 2 and 1 and reducing the remaining nodes by 2. 
This transformation is obviously a bijection. Thus, %2n+i,2,k = ^2n-i,,,k-2 
for 3 < /c < 2n — 1. The second row of Mat(T2n+i) is then equal to the 
sequence 

(5.1) 0, 0, T2n-l,.,l7 ^2n-l,,,2i ■ ■ ■ 5 '?2n-l,.,2n-37 2^271-1, ., 2n-2 (= 0); 

which is also equal to 

(5.2) 0, T2n-l,l,.(= 0), 'X2n-1,2,., ■ ■ • , 2^2n-l,2n-2,., 0, 

by Poupard's result (1.5) (also by our combinatorial proof sketched in the 
Introduction). 

(/4) The set T2n+i,m,i is empty when m = 1, 2 or 2n. When 3 < to < 
2n — 1 the diagram 



2n+l , , , , 

^m — 1 




1 



serves to illustrate the transformation that maps each tree from %2n+i,m,i 
onto a tree from %2n-i,m-i,,, by deleting the two nodes (2n-|-l) and 1, 
and reducing the remaining nodes by 1. Thus, the first column of M„ is 
equal to sequence (5.2), when read from top to bottom. 
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For the second column we first note that '^2n+i,m,2 is empty when 
m = 1, 2 and n > 2. When m = 3, the mapping 





2n+l 
1 

shows that ^2^+1, 3, 2 is in bijection with l2n-i,,,i = '2^2n-i,2,.- When 
m = 4, the following decomposition prevails 




2^271+1,4,2 — 

the two sets on the right-hand side being in bijection, respectively, with 




that is, 

= 2^2n-l,3,. + 2^2n-l,2,.- 
To prove the identity l2n+l,m,2 = '^2n-l,m-l,, + '^2n-l,m-2,, = 

'^2n+i,m,i + '^2n+i,m-i,i (by (5.2) and (/4)) for 5 < m < 2n proceed 
by induction on m using relation (-R1) already proved in (3.3). Write 

2^2n+l,m,2 = 2 T2n+l,m-l,2 — 2^2n+l,m-2,2 " 2 T2n-l,m-2,2 [by (3.3)] 

2 l2n+l,m-l,2 = 2 l2n+l,m-l,l + 2 l2n+l,m-2,l [by iuduCtioU OU m] 

- 2^2n+l,m-2,2 = -'^2n+l,m-2,l - '^2n+l,m-3,l [by iuductioU OU m] 

~ 2 'X2n-l,m-2,2 = ~2 T2n-l,m-2,l "" 2 %2n-l,m-3,l 

[by induction on n and m] 

2 2^2n+l,m-l,l — 2^271+1, m-2,1 — 2 2^271-1,771-2,1 = "^277+1,771,1 [by (3.3)] 

2 2^277+1,777-2,1 —2^277 + 1,777-3,1 ^ ^ %2n-l,Tn-3,l = 2^277+1,777-1,1- [by (3.3)] 

By summing those six equations we get 2^277+1,777,2 = 2^277+1,777,1 + 
2^277+1,777-1,1, also equal to 2:277-1,777-1,. + 2:277-1,777-2,. by (5.2). □ 

Remark. It would be interesting to make up a proof of Property 5.1 
that would have no recourse to a recurrence argument for to > 5 as above. 
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6. Proofs of Theorems 1.1 and 1.2 

Taking Property 2.1, Corollaries 3.2 and 4.2, Property 5.1 into account 
we conclude that the sequence of matrices Mat(T2n+i) is both a Delta and 
Gamma sequence. Those sequences are then identical and we may write 

(6.1) fnim, k) = %2n+l,m,k 

for all n, m, k. This completes the proof of Theorem 1.1. 

We now exploit the properties of the strictly ordered binary trees to 
prove that the matrices M^ are symmetric with respect to their counter- 
diagonals (Theorem 1.2). First, the symmetry property is banal for Mi, 
M2. For n > 3 consider the NE- and SW-corners 

/^(l,2n-l)/^(l,2n)\^ / 

/n(2, 2n - 1) U{2, n) J \U2, 2n - 1) 

7n(2n - 1, 1) Ui2n - 1, 2)\ ^ f U{2n - 1, 1) /,(2n - 1, 2) 
/^(2n,l) /n(2n,2) J \ /n(2?i,2) 

of the matrix M„. As/n(2n-l,l) = /^(2n, 2) = T2n-3 = Tsn-s/S"-^ (by 
combining (1.11), (/2), (/3), (5.1) and (5.2)), both corners are symmetric 
with respect to their counter-diagonals [in short, counter-symmetric]. 

Let us prove that the upper part of the matrix M^ is counter-symmetric 
and for z = 1, 2, . . . , n — 1 adopt the notation: 

Rowi = {(z, i + 1), (z, z + 2), . . . , (z, 2n - i), (z, 2n - z + 1)}; 
Col2n-*+i = {{i, 2?i - z + 1), (z + 1, 2n - z + 1), 

. . . , (2?i - z - 1, 2n - z + 1), (2?i - z, 2n - z + 1)}. 

Note that Row^ and Col2n-i+i have the cell (i, 2n — i + 1), belonging to 
the counter-diagonal, in common. There is nothing to prove for the pairs 
(to, k) along the counter-diagonal and also for the entries from Rowi and 
Col2n5 which are all zero. 

Let (too,/co) belong to Rowj for some j such that 2 < j < n — 1. 
Further, assume that (1.12) holds for all (m, k) G Rowi U ■ ■ ■ U Rowj_i 
and all (to, k) G Rowj lying on the right of (toq, /cq), not including (toq, /cq), 
that is, TO = j and k > kg. By symmetry, (1.12) also holds for all 
(?7i, k) G Col2n U ■ ■ ■ U Col2n-j+2 and all (m, k) G Col2n-j+i lying above 
(2n + 1 — koj2n + 1 — mo) not including the latter pair. 

Now, the following relations hold: 

/n(mo,/co) = 2/n(TOo,/co + 1) - fn{mo,ko + 2) - 2 /n-i(mo, /co) [by {R2)] 
/„_i(mo, ko) = fn-i{2n - 1 - /cq, 2n - 1 - toq), [by induction on n] 

/^(too, ko + 1) = fn{2n - ko,2n + l- toq), 
fn{mo,ko + 2) = /Tj,(2n — 1 — /co, 2n + l — mo).[both by the local induction] 
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Therefore, 

/n(mo, ko) - 2 fn{2n -ko,2n + l- mo) 

+ fnC^n — 1 — /cq, 2n + 1 — Too) + 2 /„_i(2n — 1 — ko,2n — 1 — toq) = 0. 

But by (-R3) written at (to, k) = {2n — 1 — /cq, 2n + 1 — toq) we have: 

fn{2n — 1 — /cq, 2n + 1 — mo) — 2 fn{2n — /cq, 2n + 1 — too) 
— fn{2n — /cq + 1, 2n + 1 — mo) + 2 /„_i(2n — 1 — /cq, 2n — 1 — mo) = 0. 

By comparing the last two equations we conclude that 

/n(mo, ko) = fn{2n - ko + l,2n + l- too), 

which means that (1.12) now holds for (to, k) = (mo, ko)- 

For the entries of M^ lying fee/ow the diagonal we proceed in the same 
manner and adopt the notation: 

Row2n+i-i = {(2n + 1 - z, z), . . . , {2n + 1 - z, 2n - i)}; 
Coli = {{i + l,i),{i + 2,i),...,{2n + l- i, i)}; 

for z = 1, 2, . . . , n — 1. Again, Row2n+i-i and Col^ have the cell (2n+l— z, i) 
in common. 

Let (too,A;o) belong to Colj for some j such that 1 < j < n — 1. 
Further, assume that (1.12) holds for all (m, k) G Coli U ■ ■ ■ U Colj_i 
and all {m,k) G Colj lying below (too, /co): not including (mo,/co), that 
is, TO > Too and k = ko- By symmetry, (1.12) also holds for all (to, k) G 
Row2n U ■ ■ ■ U Row2n-j-i-2 and all (to, k) G Row2n+i-j lying to the left of 
{2n + 1 — ko,2n + 1 — toq) not including the latter pair. 

Now, the following relations hold: 

/„(too, ko)=2fn{mo + 1, ko) - fn{rno + 2, ko) - 2 /„_i(mo, ko), 

[by(i?l)], 
fn-i{mo,ko) = fn-i{2n- l-ko,2n-l -mo), 

[by induction on n], 
fnimo + 1, ko) = /„(2n + 1 - /cq, 2n - too), 
fnimo + 2, ko) = fn{2n + 1 - /cq, 2n - 1 - mo), 

[by the local induction]. 
Therefore, 

fnimo, ko) = 2 fn{2n + 1- ko,2n- mo) 
- fn{2n + 1 - ko,2n-l - mo) - 2 /„_i(2n - 1 - /co, 2n - 1 - too). 
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But by (-R4) written at (to, k) = {2n + 1 — /cq, 2?i — 1 — mo) we have: 
fnC^n + 1 — ko,2n — 1 — mo) = 2 fn{2n + 1 — ko,2n — mo) 

— fn{2n + 1 — /co, 2n + 1 — mo) — 2 fn-i{2n — 1 — /co, 2?i — 1 — mo). 
By comparing the last two equations we conclude that 

/2n+i(mo, ko) = /2n+i(2?i + 1 - /cq, 2n + 1 - mo)), 
which means that (1.12) now holds for (m, k) = (mo, ko). [] 

Define 

(6.2) Eoc(t) := eoc(t), but Pom(t) := 2?i + 1 - pom(t). 
Theorem 6.2. Let 'Z2n+iix,y) := ^ ^Eoc(t)^Pom(t) 1)q tiie generating 
polynomial for the set '^2n+i by the pair of statistics (Eoc, Pom). Then, 

(6.3) %2n+l{x,y) ^%2n+l{y,x). 

Proof. This is a simple consequence of Theorem 1.2: let gn{m,k) := 
#{Eoc = m. Pom = k}. Then, 

(7^(m, k) = #{eoc = m, pom = 2n + 1 — k} = fnirn, 2n + 1 — k) 

= fn{k,2n + l-m) [by (1.12)] 

= gn{k,m). 
Thus, 

l2n+i{x,y)= Yl x'^^^'^y^^^^'^ = Y,9n{m,k)x^y^ 

t6l2n+i m,k 

^J29nik,m)x'^y'' = l2n+i{y,x). D 

Tn,k 

7. Further properties 

Other properties of the Delta Sequence can be obtained by having a 
further look at the geometry of the strictly ordered binary trees. The sub- 
and superdiagonals of the matrices M^ for n = 2, 3, 4, 5 are equal, as can 
be seen in Fig. 1.2. For an arbitrary n > 2 we have the following. 

Property 7.1. Sub- and super diagonals are equal: 

(7.1) fnik + 1, A;) = U{k, k + 1) (n > 2; 1 < A; < 2n - 1). 

Proof. First, note that k and (k + 1) can be siblings in a tree from 
'^2n+i,k,k+ij but never in a tree from %2n+i,k+i,k- Second, k can be parent 
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of (/c + 1) in a tree from the latter set, but never in a tree from the former 
one. Also, /n(2, 1) = /^(1,2) = for n > 2 and for k > 2 we have the 
decompositions: 

p, 2n+l Q 2n+l 

k \/ JO 



fn{k,k + l 





k]- 



U{k + l,k)= I + [ I , • (A; + l) 

The first terms in the previous two equations are in bijection, as well as the 
second ones, the notation "• k" meaning that k is the end of the minimal 
chain, following our convention on Tree Calculus (c/. Section 3). [] 

Property 7.2. We have the crossing equahties: 

(7.2) f^{k + l,k-l) + Uk-l,k + l) 

= U{k + l,k) + U{k-l,k) 

= U{k,k + l) + U{k,k-l) {2<k<2n-l). 

Reporting to Fig. 1.2 the involved entries in the first identities are 
located on the four bullets drawn in the following diagram. 

k-1 k k + 1 

k-1 

k 
k + 1 

Proof. Let i, j he two different integers from the set {(/c — l),/c, (/c+1)}. 
Say that i and j are connected in a tree t, if the tree contains the edge 
i — j, or if i and j are brothers and one of them is the end of the minimal 
chain of t. Each of the four ingredients of the previous identity is now 
decomposed into five terms, depending on whether the nodes (/c — 1), /c, 
{k + 1) are connected or not, namely: no connectedness; only k, {k + 1) 
connected; (/c — 1), /c connected; (/c — 1), (/c + 1) connected; all connected. 
Thus, 

2n+l I — I 




fn{k 




k-1 



r ,k + i] + [ 




2n+l 



k-1 




A1 + A2 + A3+A4 + A5; 
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2n+l 




B1+B2+B3 + B4+B5] 

2n+l 
fc+1 





fc-1 



.R tV° 






k+l 

:=Ci + C2 + C3 + C4 + C5; 

2n+l 




fc-l 

/„(A;-1,A; + 1) = [ 





fc-i 



:=Di + D2 + Ds + D^ + D5. 

Now, the following identities hold: Ai = Ci, A2 = C3, A3 = D4, 
A4 = Ds, 2 A5 = 1^5; Bi=Di,Bs = D2, B5 = C5. Moreover, 



fc+i 



Ba — Co 



Ca — Bo 




D^ = 2 Af 



A.. 



Altogether, E.(^^ + B,) - E.(C^ + 1^0 = (S4 - C2) - (C4 - S2) + (A5 
L>5) = 2 A5-A5 + (A5- 2^5)^0. D 
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Property 7.3. The row sums /^(m, .) form a Poupard Triangle, the 
initial conditions being: /o(l,.) = 1, /^(l,.) = and /^(2,.) = 
2 ^^ /^_i(m, .) (n > 1); and the Gnite difference system: 

(7.3) A2/n(m,.)+2/„-i(m,.)=0 (1 < m < 2n - 1). 

The column sums fni'^k) form a Poupard Triangle, the initial conditions 
being: /o(.,0) = 1, /„(.,0) = and /„(.,!) = 2 Efc/n-i(.,fc) (n > 1); 
and the finite difference system: 

(7.4) A2/,(.,fc)+2/„_i(.,A;)=0 (0<A;<2n-2). 

k 

There are several proofs of this Property. First, the methods developed 
in Section 3 can be readapted by disregarding the conditions involving the 
pom-statistic. Here, we simply work out a specialization of the recurrence 
relations (-R1) — (-R4), that makes use of the previous two properties. 
Besides, we only prove the first part of the property that deals with the 
row sums. 

Proof. For 1 < m < 2??, — 2 we have: 
A2/n(m,.)= V A'fn{m,k) 

m ' ' m 

l<k<2n 



V A'Um,k) 

< ' m 

t<m — 1 

+ A^{fn{m,m) + fn{m,m + l) + /n(m, m + 2)) 



m 



+ V AVn(m,fc) 

' ' m 

m+3<k<2n 

= ^ -2/^_i(m, /c) 

l<fc<m-l 

+ /n(m, m) - 2fn{m + 1, m) + fn{m + 2, m) 
+ fn{rn, m + 1) -2fn{m + l,m + l) + fn{m + 2,m + 1) 
+ fn{m, m + 2) - 2/^(m + 1, m + 2) + /^(m + 2, m + 2) 
+ Yl -2/n-i(m,A;-2). 

m+3<k<2n 

In the previous sum the diagonal terms vanish. Also, fni'tn + l,m) = 
fn{m,m, + 1), fni'tTi + 2, ?7i + 1) = fni'tTi + l,m + 2) by (7.1). The sum 
of the nine intermediate terms becomes: —fnijn + 1, m) + fnijn + 2, m) — 
fn{m + l,m + 2)+ fn{rn,m, + 2), which is by (7.2). Hence, 

A2/n(m,.) = -2 V U-i{m,k)-2 V U-i{m,k) 

l<fe<m-l m+l<k<2n-2 

= -2/„_i(m, .). 
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On the other hand, by (/2) and (1.11), 

/„(2n-l,.)= Yl fn{2n-l,k)= Yl (/n-i(A;, .) + /n-i(., A;)) 

l<k<2n l<k<2n-2 

= 2 /n-l(», •) := 2 #l2n-l = 2T2n-l/2"' ; 

/^(2n,.)= Yl f'-{2n,k)= Y /-i(^'-) 

l<k<2n l<k<2n-2 

= /n-l(», •) = #2^2n-l = 72^-1/2" ; 

/„(2n + l,.)=0. 



The last three evaluations imply that /^fn{'2.n — 1, 

m 

i,.)-o. D 



-2/„_i(2n 



8. Other equivalent definitions for the Delta sequence 

Definitions 1.1 and 1.2 have been shown to be equivalent to characterize 
a Delta Sequence. Other combinations of the recurrence relations (-R1)- 
(i?4), together with the initial conditions (/l)-(/4), can be used. We 
describe them by means of squares and arrows, as was done in Fig. 1.3: 



(i?3) 



RoAVi 
R0W2 



Row2„-i 

R0W2n 





- 


\ 






\ 




\ 



C0I1C0I2 Col2„_l 

rq nCoi2„ 



(i?4) 4- 




>(i?2) 



(i?l) 
Definition 1.3 



Definition 1.4 



Other initial conditions than (/I) — (/4) can be introduced. They will 
be denoted by {SWl.A), {NE2,3), as they refer only to the South- West 
and North-East corners of the matrices: 



{SW) 
(NE) 



/„(2n-l,l) /„(2n-l,2) 
/„(2n,l) /„(2n,2) 



/n-i(.,l) /„-l(2,.) + /n-l(.,2) 
/™-i(.,l) 



/„(l,2n-l) Ml,2n)\^( 0\^ / 

/„(2,2n-l) /„(2,n) ; \U-ii2,,)0j {U-ii^n - 2, ,) 



When one of those two conditions (SW), (NE) is involved, two recur- 
rence relations among {R1)-{R4) are needed to build up an equivalent 
definition. In Definition 1.5 for instance, (-R1) and (-R4) are to be associ- 
ated with (SW). We then get five further equivalent definitions: 
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(i?3) 



(fl3) 




Col2„_l 

Col2„ 




>(i?2) 



Row2„_i 
Row2„ 




(Rl) 
Definition 1.5 



(Rl) 
Definition 1.6 



(Rl) 
Definition 1.7 



(i?3) 
Coll Cob 



(i?3) 



(i?4) i 




Definition 1. 



(Rl) 
Definition 1.9 



We do not reproduce any proofs for those equivalences, but point out 
the fact that our Tree Calculus requires that each initial condition be 
combinatorially interpreted, as was done in Sections 2 and 5. 

9. Generating functions for the Delta sequence 

1. Poupard matrices. Let G = (gij) {i > 0, j > 0) be an infinite 
matrix with nonnegative integral entries. Say that G is a Poupard matrix, 
if for every i > 0, j > the following identity holds: 



(9.1) 



9i,j+2 — 2 Qi+ij+i + gi+2,j + 2 gij — 0. 



Let G{x,y) := E 9r,j (^7^) (?/Vj!); Mv) ■= E 9r,j (j/VjO (^ > 0); 
j>o,i>o j>o 

Cj{x) := Yl 9i,j (x'^/i) {j > 0) be the exponential generating functions for 

i>0 

the matrix itself, its rows and columns, respectively. 

Proposition 9.1. The following four properties are equivalent. 
(i) G = (gij) {i > 0, J > 0) is a Poupard matrix; 
(ii) R'l{y) - 2 R'^^M + R^+2{y) + 2 i?,(y) = for all i > 0; 
(Hi) C'^{x) - 2 q+i(x) + C,+2(x) + 2 Cj{x) for all j > 0; 
(iv) G{x,y) satisfies the partial differential equation: 



dx'^ 



dxdy 



dy"' 



Proof. It suffices to write R'l{y) - 2R[^^{y) + Ri+2{y) + 2i?i(y) = 
E(fi'i,i+2 - 'igi+id+i + 2fi'i+2,i - 2g,j){y^/j\), and, similarly, Cj{x) - 

j>0 
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2Cj^^{x) + Cj+2{x) + 2Cj{x) = Y.{9i+2,j - 2^i+ij+i + 2^ij+2 - 

i>0 

2gij){x'^/il) to obtain the equivalence between the first three properties. 
As for the last one, simply note that G{x,y) = X]i>o -^^(2^) -^V^' = 
Y^ Cj{x)y^ /jl and make the appropriate derivations. [] 

j>0 

Proposition 9.2. We have 

(9.3) G{x, y) = A{x + y) cos{V2y) + B{x + y) sin(V2 y), 

where A{x) and B{x) are two arbitrary series. 

Proof. Let ^ := x + y, rj := y. Then, dG/dx = dG/d^; dG/dy = 
dG/dC + dG/drj; d^G/dx^ = d^G/d^; d^G/dy^ = d^G/d^^ + 
2 d^G/dC drj + d^G/drj^; d^G/dx dy = d^G/d^^ + d^G/di drj. 

Thus, (9.3) can be rewritten as 

(9.4) + 2G = O. 

The ordinary differential equation G" + 2G = 0, whose characteristic 
polynomial is r^ + 2 = 0, has a general solution of the form A cos{\/2r)) + 
B sin(v2?7), so that the general solution of (9.2) is exactly given by 
(9.3). D 

The exact expression for the generating function G{x,y) can then be 
derived, if A(x + y) and B{x + y) can be obtained by an independent 
calculation, as done in the sequel. 

2. A sequence of Poupard matrices for the lower triangles. The entries 
/n(m, k) {1 < k < m < 2n) from the lower triangles in the matrices M^ 
(n > 1) are now recorded as entries A^^- (p > 0, z > 0, j > 0) of infinite 

matrices A^^^ = (A^^^) (z > 0, j > 0) as follows. 
Define 



eg 5N ^{p) .= / 0. ifi+j=p mod 2; 

^ ' ' ^'^ ' |^/n(m, /c), if z + J = p + 1 mod 2; 

with k := j'+l, m := i+j+2, 2n := p+{i+j)+l. The latter equation makes 
sense, as i+j and p are of different parity. The mapping (p, z, j) \-^ (n, to, k) 
is one-to-one, the reverse mapping being: p := 2n — m + l, i := m — k — 1, 
j :— k — 1. Thus, for p > 0, the matrix A'-^^"'"^-' reads: 
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//p+i(2,l) 



1 
2 
3 

4 

l\ 



/p+2(4,3) 



/p+2(4,2) /p+3(6,4) 

/p+2(4,l) /p+3(6,3) 

/p+3(6,2) 

/p+3(6,l) 





4 5 

/p+3(6,5) 



and for p > 1, the matrix A^'^^' is equal to 
12 3 

/ /p+i(3,2) /p+2(5,4) 



1 

2 
3 

4 



/p+i(3,l) /p+2(5,3) /p+3(7,5) 

/p+2(5,2) /p+3(7,4) 

/p+2(5,l) /p+3(7,3) 

/p+3(7,2) 

/p+3(7,l) 



/p+3(7,6)\ 



/ 



Proposition 9.3. Every matrix A^P' {p > 0) is a Poupard matrix. 
Proof. Using Definition (9.5) we have 

= /n+i(m + 2, A; + 2) - 2 U+i{m + 2, A; + 1) 

+ /„+i(m + 2,A;) + 2/„(m,/c) 
= A /n+i(m + 2, /e) + 2 /„(m, A;) = 0, 

k 

by rule (i? 4). □ 

3. The first matrix A^^\ The counter-diagonal of A*^^^ at depth 2i 
(z > 0) reads 



AlP 



.(1) 



2, n - /^+i(2z + 2, 1); A^^i, 1 = /,+i(2z + 2, 2); 



An,- 



(1) 



2J-j,i 



/,+i(2z + 2,j + l) 



a[,;^, = /,+i(2z + 2,2z + 1) 



Those (2i-(-l) terms are equal to the first (2z + l) entries of the (2i-|-2)-nd 
row of the matrix M^_|_i, that is, to the single term /i(2, 1) = 1 for i = 
and for i > 1 to 



/,(1,.)=0, /.(2 



5 '^5 • • • 5 



Mj + i 



7 ';? • • • 5 



/^(2z,.),0, 



by virtue of (/2). Thus, A*^^) is identical to the Poupard matrix (/^(j+l, .)) 
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equal to 








1 


2 


3 


4 5 


6 





(l 











o\ 


1 


/i(2,.) 





/2(4,.) 


/3(6, 


.) 


2 





/2(3,.) 





/3(5,.) 




3 


/2(2,.) 





/3(4,.) 






4 





/3(3,.) 








5 


/3(2,.) 










6 


\o 








/ 





1 2 


3 4 


5 6 ••• 







( 1 








--A 




1 





1 


1 


4 ... 




2 





2 


8 






3 





1 


10 ■•• 






4 





8 






) 


5 





4 ... 










6 














4. The first two columns of A^p\ In the sequel the column labeled j 
of A^P^ will be denoted by A^^- and the exponential generating func- 
tion for that column by A]^Ux) = X]i>o -^i i -^V^'- Also, A^P\x,y) := 
^jyQAl^Ux)y^/j\ will be the double exponential generating function for 
the matrix A^^'. 

Proposition 9.4. The first two columns of each, matrix A^p^ {p > 2) are 
related to the columns of A^^^ by the identities: 

(9.6) Ai^)(a:) = Ai]l_,{x), Ai^J(x) = ^^^^,l-i{x) + A^^^lix). 

Proof. For the first identity it suffices to prove A^q = A^ _j^, that is 

/,(z + 2,l) = /,(p + z + l,p) 

when 2n = p + i + 1. This is true by Theorem 1.2. For the second identity 
it suffices to prove X^^ = A^_^\ p-i + -^i p? that is 

(9.7) fnii + 3, 2) = fn{p + i + 2,p) + f^{p + i + 2,p+l) 
when 2n = p + i + 2. But by (/ 4) 

/„(z + 3,2) = ^(z + 2,l) + /„(z + 3,l), 
so that by Theorem 1.2, identity (9.7) holds. [] 
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5. Solving the partial differential equation. By (9.6) 



^/'\x,y) 



{y=o} 



{y=o} 



Ai"(l(a:) = Ai;i_i(x); 



As each hS'^' is a Poupard matrix, we can use identity (9.3), so that 



A(p)(x,2/) 



{y=o} 



A{x) and dA^P\x,y)/dy = A'{x + y) cos{y/2y) - 



y2A(a; + 2/) s\n{^/2y)+B'{x + y) siiv{^/2y) + ^/2B{x+y) cos (v^y). Hence, 



A(^)(x,y) 

M(p)(x,y) 
dy 



{y=o} 



{y=o} 



^(^) = A;!^_i(x); 



A'(x) + V2S( 



X) 



d (1) 



■Aili(x) + AW(a:). 



Consequently, A(x) = h.\ ^^{x) and i?(x) = K,^p{x)/\/2 and the general 
expression for KP{x,y) reads: 

(9.8) A(p)(x, y) = A[]1_^{x + y) cos{V2y) + AW(x + y) sin{V2y)/V2. 

This expression still holds for p = 1. We know that A^ q(x) = 1. On 
the other hand, the coefficient of x'^^^^ /{2k + 1)! {k > 0) in A^ {{x) is 
fk+i{2, .) = T2k+i/2K Hence, A;\^(x) = ^2 tan(x/y2) and 



A(^Hx,y) 



cos 



(9.9) 



cos 



X — y\ I (X + y 



V2 



cos 



v^ 



a result already obtained by Poupard. 

Remark. For getting the solution for A*^^''(a:, y) we can also start with 
the general expression displayed in (9.3) and calculate A and B with the 
initial conditions A(^)(0,y) = A(^)(a;,0) = 1. We find 



(9.10) 



A(^)(x,y) = cos(y2,) + ^--^<f(^ + y)) 3i,(^^) 

sin(v2(x + y)) 

sin(v^a^) + sm.{\/2y) 



sin 



{V2{x + y)) 



an expression, which is naturally equal to the right-hand side of (9.9) (by 
a simple trigonometric calculation) . 
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We have not worked out other expUcit formulas for AS'P'>{x^y) when 
p > 3, but only derived the exponential generating function for those 
series, as explained in the next subsection. 

5. A generating function for the lower triangles. By (9.9) 
Iv-^'{x,y) = ^:^ cos( — ^) /cos 



dy ' '"' dy V ^2 / / V ^2 

'X + y 

V2 



(9.11) = sin(y2x) / V2 cos^ (- 

LetA(x,y,z):=5^A(p)(x,y)^^^. 



Then, 



^"-' .(1) 



A(x, y, z) = cos(V2 y) J^ (j, _ i)i Kp-i(^ + v) 

p>i ^^ '' 



+ "^^4^ > : zf— rT7Ai:^(^ + 1/) [by (9.8)] 






cos(v^y) A(i)(x + y, .) + ^^^^^^A(i)(x + y, .), 



since K^^\x,z) = ^ — Ai^^(a;). 
p>0 ^' 
By (9.11) we then get: 

A(x,y,z) =cos(v^y) cos(^ -= — j / cos(^ -= — 



2fx + y + z 



+ sin(y2y)sin(y2(x + 2/)) /2 cos^, 

' ^ v2 

cos(-\/2x) + cos(v^y) cos(v^ 2) 

~ 2fx + y + z 
2 cos^ 



V2 
Now, express A(x, y, 2) as a series in the /n(?Ti, /c)'s. By definition, 

A(x, ,, .) = J2 a15^^|| (p > 1, z > 0, , > 0); 
so that by (9.5) 

™m — fc — 1 „ k — 1 y2n—m 

A{x,y,z)= V fn{m,k)- ^ — — — — — -, 

^^-^ [m — k — 1); [k — Ij! (zn — m)! 
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the latter sum over the set {2<A; + l<m< 2n}. This achieves the proof 
of Theorem 1.4. 

Note that the right-hand side of (1.14) is symmetric in y,z, that is, 
A{xjy,z) = A(x,z,y). The change y ^ z in the left-hand side of (1.14) 
shows that 

fn{2n + 1 - k,2n + 1 - m) = fn{rn, k), 

the symmetry proved for the entries /n(?Ti, k) such that m > k + 1. 
Remark. Let ^ = in (1.14). We get 



2n-k-l yk-l C0s(y2x)+C0s(y2y) 



:+l<2n ^ J \ J 2 COS^ ■:=- 

V ./2 



or 

\^ n(i)^'^' _ ^{l)t^ ..^_ cos(v^x) + cos(v/2y) ^ 

i>o,i>o -^ 2 cos^l — ^1 

which is another expression for A^^^x^y) than (9.9) and (9.10). 

6. A generating function for the upper triangles. The entries fnitn, k) 
{1 < m < k < 2n) from the upper triangles in the matrices M^ (n > 1) 
are next recorded as entries a;^ (p > 0, z > 0, j > 0) of infinite matrices 
0(p) = (a;g)) (i > 0, J > 0) as follows. 
Define 

jp) .^ f 0' if « + J ^ P mod 2; 

''^ ' Xfnirn.k), ifz+j=pmod2; 

with m := p + 1, k := p + j + 2, 2n := p + i+j + 2. Conversely, i := 2n — k, 
j := k — m — 1, p := m — 1. Thus, for p > 

Q(2p+1) 



1 2 3 





/ /p+2(2p + 2,2p + 4) •••\ 

1 /p+2(2p + 2,2p + 3) /p+4(2p + 2,2p + 5) •• 

2 /p+3(2p + 2,2p + 4) 

3 /p+3(2p + 2,2p + 3) 
4V ••• / 

^^ 1 2 3 

//p+i(2p+l,2p + 2) /p+2(2p+l,2p + 4) o\ 

1 /p+2(2p+l,2p + 3) 

2 /p+2(2p+l,2p + 2) 

3 
4\ ■■■ / 
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Remark. The first rows of all matrices M„ are null, so that O*-*^^ is the 
infinite matrix with all entries equal to zero! 



Also 


, write 














n(i) = 












1 


2 


3 


4 


5 6 7 





( /2(2,4) 





/3(2,6) 





/4(2,8) /5(2,10)\ 


1 


/2(2,3) 


/3(2,5) 





/4(2,7) 


/5(^ 


^,9) 


2 


/3(2,4) 





/4(2,6) 





/5(2,8) 




3 


/3(2,3) 


/4(2,5) 





/5(2,7) 






4 


/4(2,4) 





/5(2,6) 








5 


/4(2,3) 


/5(2,5) 










6 


/5(2,4) 












7 


\/5(2,3) 










/ 


12 


3 4 


5 6 


7 









/ 








o\ 








1 


1 1 


4 


34 










2 


2 


8 


68 ••• 










3 


1 10 


94 












4 


8 


104 • • • 












5 


4 94 














6 


68 • • • 
















7 


^34 ••• 






/ 









Proposition 9.6. Every matrix 0*^^-' (p > 0) is a Poupard matrix. 
Same proof as for Proposition 9.3. 
The row labeled i of fl^^' will be denoted by O^^J and the exponential 



(p), 



,(p). 



generating function for that row by fl\^{y) = ^j>o^i Jy-' /j^-- Also, 
Q^P\x,y) := ^i>o^i,{y)x^/i^- will be the double exponential generating 
function for the matrix fl^^'. As x and y play a symmetric role in (9.2), 
the solution in (9.2) may also be written 

G{x, y) = A{x + y) cos(v^x) + B{x + y) sin(v^ x), 
so that the generating function of each matrix VL^P' is of the form 

Vt^p\x,y) = A{x + y) cos(V2x) + B{x + y) sin(y2a:). 

As the first row of each matrix Q,^p' is the zero sequence, we have 
O(p)(0,y) = A{y) = 0. Hence, n'^P^x.y) = B{x + y) sin(y2x) and 



d 



d 



n^P'^ (x, y) = (-K-B{x + y)) sin( V2x) + V2 B{x + y) cos(y2; 



X) 
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Therefore, 



dx 
and then 



{x=0} 



(9.12) n^P\x, y) = ^ sin(V2x) fl'flix + y). 

v2 



The evaluation of Q\ l(y) is easy, as the row labeled 1 of Q^^' is 
(1, 0, 1, 0, 4, 0, 34, 0, ... ), compared with (0, 1, 0, 1, 0, 4, 0, 34, 0, ... ), which 
is the sequence of the coefficients of the Taylor expansion of y^ tan(y/v^). 
In fact we have u[]lj = /i+2(2, 2j + 3) = fj+i (2j + 2, .) = 72^+1/2^. Thus, 

d 



"•-<")- 1; ^-B 



dy \^J cos2(y/v^) 

so that 



(9.13) 0(^)(x,2/) = ^sin(V2 



1 

X) 



a/2 2{ ^ + y \' 

V COS I — -j=- I 



Proposition 9.7. For all p > 1 we have: 

(9.14) <^(y) = oW(y)- 

Same proof as for Proposition 9.4. Now, define: 
(9.15) 0(x,2/,^):=5^0(p)(a:,2/)^ 

and make use of (9.18)— (9.20): 

1 P 

n{x,y,z) = -r=sm{V2x)^n[^^,{x + y) — 

^l=sm{V2x)J2^l]l{x + y)^ 
^^ p>i ^' 



^ sin( V2 x) O^^^ (2, a; + y) 
sin(v2a;) sin(v22) — 



2/^^ + 2/ + ^ 

2 cos ' 



V2 
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TREE CALCULUS FOR DIFFERENCE EQUATIONS 
As all the entries UqJ (p > 1, J > 0) are null, 

n{x,y,z) = E'^SStT^T iP>h^>hJ> 0); 

P,r,J 

(p) 

SO that by definition of the ooljS 

™2n — A: „ k—m — 1 ^m — 1 

(9.16) n{x,y,z)= V fn{m,k)— — — — — , 

■^^ [In — k)\ [k — m — 1)\ [m — Ij! 

the latter sum over the set {2 < m + 1 < k < 2n — 1}. This achieves the 
proof of Theorem 1.5. 

The right-hand side of (1.15) is symmetric in x, z, that is, Q{x, y, z) = 
Q{z,y,x). The change x ^ z in the left-hand side of (1.15) shows that 

fn{2n + 1 - k,2n + 1 - m) = fn{rn, k), 

the symmetry proved for the entries fniin, k) such that m + 1 < k. 
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